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Abstract 

In the context of pure QED, we obtain analytic expressions for the contributions 
to the Bhabha scattering differential cross section at order which originate from 
the interference of two-loop photonic vertices with tree-level diagrams and from the 
interference of one- loop photonic diagrams amongst themselves. The ultraviolet 
renormalization is carried out. The IR-divergent soft-photon emission corrections 
are evaluated and added to the virtual cross section. The cross section obtained 
in this manner is valid for on-shell electrons and positrons of finite mass, and for 
arbitrary values of the center of mass energy and momentum transfer. We pro- 
vide the expansion of our results in powers of the electron mass, and we compare 
them with the corresponding expansion of the complete order photonic cross 
section, recently obtained in ^U]. As a by-product, we obtain the contribution to 
the Bhabha scattering differential cross section of the interference of the two-loop 
photonic boxes with the tree-level diagrams, up to terms suppressed by positive 
powers of the electron mass. We evaluate numerically the various contributions to 
the cross section, paying particular attention to the comparison between exact and 
expanded results. 
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1 Introduction 



The Bhabha scattering process plays a crucial role in the study of elementary particle 
phenomenology, since it is the process employed in the luminosity measurement at e~^e~ 
colliders. The small-angle Bhabha scattering at high (~ 100 GeV) center of mass energy 
and the large-angle Bhabha scattering at intermediate (1 — 10 GeV) energies have cross 
sections that are large and QED dominated; these two characteristics allow for precise 
experimental measurements, as well as for a detailed theoretical evaluation of the cross 
section. 

The radiative corrections to the Bhabha scattering in pure QED have been extensively 
studied (see and references therein). The 0{a^) corrections have been known for a 
long time, even in the full Electroweak Standard Model [2]. The second order QED 
corrections, C(a'^), were the subject of renewed interest in the last few years, and several 
works have been devoted to the study of second order radiative corrections, both virtual 
and real, enhanced by factors of ln'^(s/m^) (with n = 1, 2, s the c. m. energy, and m the 
mass of the electron) |S1 IH E] • The complete set of these corrections was finally obtained 
in ini by employing the QED virtual corrections for massless electron and positrons of 
[7j, the results of jH], and by using the known structure of the IR poles in dimensional 
regularization Very recently, the complete set of photonic 0{a'^) corrections to the 
cross section that are not suppressed by positive powers of the ratio m^/s were obtained in 
|10j . The subset of virtual corrections of 0{a^) involving a closed fermion loop, together 
with the corresponding soft-photon emission corrections, were obtained in an analytic 
and non approximated form in ^T] . The results contained in these papers do not 

rely on any mass expansion, and they are valid for arbitrary values of the c. m. energy 
s, momentum transfer t, and of the electron mass m. In ^ll| IT^ IT!^ . the calculation of 
the relevant loop diagrams was performed by employing the Laporta-Remiddi algorithm 
[Ti] which takes advantage of the integration by parts ^2] and Lorentz-invariance [TBj 
identities in order to reduce the problem to the calculation of a small set of master 
integrals. The master integrals are calculated using the differential equations method 
[TYj : their expression is given in terms of harmonic polylogarithms ^H]- Both IR and 
UV divergencies were regularized in the dimensional regularization scheme JH] and they 
appear, in the intermediate results, as singularities in {D — 4), where D is the dimension of 
the space-time. It is natural to apply the same approach to the calculation of the complete 
set of 0(a^) virtual corrections. At present, the list of the master integrals required to 
complete the calculation is available in [201 1^ • However, only very few of the master 
integrals related to the two-loop photonic box diagrams have thus far been calculated 
[221 1201 1^- The subset of second order radiative corrections due to the interference of 
one-loop diagrams was studied in [21]. 

In using the results of [^1221; we calculate in this paper the following O(a^) correc- 
tions to the differential cross section in QED: 

• corrections due to the interference of the two-loop photonic vertex diagrams with 
the tree-level amplitude; 

• corrections due to the interference of one-loop photonic diagrams amongst them- 
selves; 



1 



• corrections due to the emission of two real soft photons from a tree-level diagram 
and of one soft photon from a one-loop photonic diagram. 

All of the contributions mentioned above are calculated for finite electron mass and for 
arbitrary values of the c. m. energy s and momentum transfer t. We employ dimensional 
regularization in order to regularize both UV and IR divergencies. The UV renormaliza- 
tion is carried out in the on-shell scheme. By following the same technique of fH], we 
pair virtual and soft-photon emission corrections in order to check the cancellation of the 
IR singularities. 

We also expand our results, that retain the full dependence on the electron mass, in the 
limit in which the electron mass is negligible with respect to the Mandelstam invariants. 
In this way, it is possible: 

• to prove that our results reproduce the correct small-angle Bhabha scattering cross 
section at 0{a'^), which is determined by the Dirac vertex form factor only 

• to provide strong cross-checks of a large part of the result of |10j : 

• to find, by subtracting our result from the cross section of ^Oj, the contribution 
to the cross section of the interference between two-loop photonic boxes and the 
tree-level amplitude; 

• to investigate the numerical relevance of the terms suppressed by positive powers of 
the electron mass. 

This paper is structured as follows: after a brief summary of our notation in Section |2l 
in Sectional we discuss the irreducible two-loop vertex photonic corrections, providing an 
expression for their contribution to the virtual differential cross section. In Sections |3J |31 
and El we calculate the interference of the reducible vertex diagrams with the tree-level 
amplitude, the interference of the one-loop vertex diagrams with themselves and with the 
one-loop box diagrams, respectively, and we obtain the corresponding contributions to the 
virtual differential cross section. In Section [71 we complete the analysis of the interference 
among one- loop diagrams by considering the interference of one-loop boxes. In Section |H1 
we discuss the soft-photon emission at 0{a^) and, in Section IHl we explicitly show how 
the cancellation of the IR divergencies works between virtual and soft corrections. In 
Section we analyze the expansion of our results in the limit m'^/s — > 0; we discuss 
the behavior of the cross section at small angle, compare our calculations with the results 
present in the literature, and discuss the numerical accuracy of the expansion. Section ITT] 
contains our conclusions. In Appendix ^ we collect the definition of some function 
introduced in the paper. Finally, in Appendix ^ we provide the expressions of the 
contribution of the two-loop photonic boxes to the differential cross section at order 
in the limit m^/s — > 0; all of the functions of the Mandelstam variables introduced and 
employed throughout the paper are available in electronic format in ji26j . 
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2 Kinematic, Notation and Conventions 



In this paper, we employ the notation and conventions adopted in [HI [T21 HB] , which are 
summarized in the present section. We consider the Bhabha scattering process: 



'{Pi) + e^{p2 



where pi, p2, Ps, and are the momenta of the incoming electron, incoming positron, 
outgoing electron, and outgoing positron, respectively. All of the external particles are 



on their mass-shell, i. e. p^ 



-m 



1, . . . , 4), where m is the electron mass. 



The Mandelstam invariants s, t, and u are related to the beam energy (E) and scat- 
tering angle in the center of mass frame of reference (6) by the relations 



u 



.p2 



-{Pl+P2y 
-{Pl-P^f 

-{Pl -Pif 



-4 {E' 
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-4 [E^ - m^) cos' 



(2) 
(3) 

(4) 



Moreover, the Mandelstam invariants satisfy the relation s + t + u = Am^. 

The analytic expressions for the the vertex and box form factors that we employ in 
the rest of the paper are calculated in the non physical kinematic region s < 0, and are 
then analytically continued to the physical region s > Am?. In order to express the results 
(for s < 0) in a compact form, it is convenient to introduce the dimensionless variables x, 
?/, and defined through the following relations: 



1 



X] 



-m 



-m 



X 



\/4 
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yf 
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u = — m 



z 
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\/Arrfi 
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\/Aiv? 




J —u 



< X < 1 



< y < 1 



< z < 1 



(5) 
(6) 
(7) 



The analytic continuation of all the form factors to the physical region s > Am? can be 
readily obtained by replacing x — > —x' + it, where e is an infinitesimal positive quantity 
and where 

a/s — Am^ — y/s 



X 



y/s — Am^ + ^/s ' 



(8) 



(see (12] for a more detailed discussion). 

The Bhabha scattering differential cross section, calculated by summing over the spins 
of the final state and averaged over the spins of the initial one, can be expanded in powers 
of the fine structure constant a as follows: 



dVL 



i=V,S 



a\ da^^\s,t,m'^) / a\'^ da2\s,t,m^) 

dn ^ wJ dn 



(9) 
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Figure 1: Two-loop irreducible diagrams contributing to the photonic vertex corrections of 
the electron current in the t-channel. The two diagrams, analogous to Figs (c) and (d), 
which have an insertion of a one-loop vertex and self-energy correction on the outgoing 
electron line, are not shown. 



where the superscripts V and S indicate virtual and soft-photon emission contributions, 
respectively, while the subscripts 0, 1, and 2 label the tree-level, 0{a^) and 0{a^) cor- 
rections, respectively. 

The tree-level and 0{a^) corrections are well known (their explicit expressions 

are collected, for example, in jT2l US])- The contribution of 0{a^) to the virtual cross 
section can be further split as follows: 

dVL ~ dVL ^ dVL ^ dVL ' 

(10) 

where the superscript Np = 1 indicates the UV-renormalized diagrams including a closed 
fermion loop (calculated, together with the corresponding soft radiation diagrams, in 
[T^ E])- The superscript "ph Vertices" indicates the contribution of all of the UV- 
renormalized photonic corrections which include at least one vertex diagram, while the 
superscript "ph Boxes" indicates the contribution to the cross sections of the photonic 
corrections which include box diagrams only. 

This paper is dedicated to the calculation of da^'^^ '^^^^^'^^^^ /dQ, as well as to the 
calculation of the contribution to da^'^^ /dQ deriving from the interference of one- 
loop box diagrams amongst themselves; the soft-photon emission corrections that cancel 
the residual IR poles in the quantities mentioned above is also discussed. By employing 
these results and the findings of ^U]; we also obtain the contribution of the interference 
of the two- loop photonic boxes to the cross section at order a^, up to terms of order m^/ s 
excluded. 



3 Irreducible Two-Loop Vertex Corrections 

In this section, we obtain the contribution of the irreducible photonic vertex corrections 
to the Bhabha scattering differential cross section at order a^. Four of the two-loop 
irreducible vertex graphs, which correct the electron current in the t-channel photon 
exchange contribution to the Bhabha scattering, are shown in Fig. d There are two other 
graphs contributing to the process: these are the specular image of diagrams [Il-(c) and 
n]-(d) with the fermionic arrow reversed, and their contribution to the differential cross 
section is identical to that of the diagrams in[T]-(c) and[T]-(d). The sum of the graphs[T]-(a) 
and[T]-(b), and of the graphs [T]-(c) and[T]-(d) are separately gauge independent. 
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(a) 



(b) 



(c) 



(d) 



Figure 2: Two-loop photonic vertex diagrams contributing to Bhabha scattering. 



The explicit expressions of diagrams n-(a)-n-(d) were calculated in [221 • The electron 
current which includes the two-loop photonic corrections can be written as 

2 



ay 

TT J 



pM^)^^ + -L^f {p, - p,) 
2m 



where a'^'^ = —i/2['~f'^, '-/"]. The electron spinors and the dependence of the form factors 
= 1,2) on the electron mass m are omitted in Eq. The expression of 



F^'^'^'^t) 

the contribution of the single diagrams to the UV-renormalized form factors F'f'^'^^\t) 
can be found in [22], and in jSEj- The form factors shown in Eq. (|TT|l still include IR 
singularities, which are regularized within the dimensional regularization scheme. The 
Laurent expansion of the form factors is 



{2l,ph,-2) 



(D-A) 



{2l,ph-l) 



(D-4) 



+ Ff'^''°)(t) + C>(D-4) 



^{2l,ph,-l) 



it) 



it) 



{2l,ph,0) 



(t) + C(D-4) 



(12) 



(D-A) 

where D is the dimensional regulator. 

The four photonic two-loop vertex correction diagrams contributing to Bhabha scat- 
tering are shown in Fig. |2l where the shaded circle represents the sum of the vertex graphs 
shown^ in Fig. ^ 

The interference of diagrams |21-(b) and|21-(d) with the tree-level amplitude provides 
a contribution to the differential cross section that is identical to the one of diagrams 
|2]-(a) and|21-(c), respectively. This can easily be proved by observing that diagram |21-(b) 
(j2Kd)) can be obtained from|21-(a) (j21-(c)) by applying the transformations p2 ^ —ps 
and p4 — Pi, and that these transformations leave the Mandelstam invariants s and t 
unchanged. 

The contribution of the diagram in Fig. 121- (a) to the Bhabha scattering differential 
cross section at order a"^ can be written as 

.,2 r 1 1 

(2«,Irr. Ver.)/ ,n J- ^ r(2i,Irr. Ver.) , 



^^(y,ph Vertices) 



(phlrr. Ver.,2-a) 



a 
s 



St ^ 



s,t) + :^V,'- 



{s,t) 



(13) 



The functions V^^'^^'^^^'^^^'\s,t) {i = 1,2) have the Laurent expansion 



(2«,Irr. Ver.) 



S,t) 



(2«,Irr. Ver.,-2) 



S,t) 



{D - 4)2 

(2/,Irr. Ver.,0) 



+ 



V^''''-''-^''\s,t) 
{D-A) 



s,t) + 0{D-4) 



(14) 



The diagrams in Fig.UKc) andEKd) enter the sum with a multiphcity factor of 2. 
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with 



{2/,Irr. Ver.,-2)j 
(2Z,Irr. Ver.,-1) , 



S,t) 



(15) 



Ca{s,t)ReFi''''''''\t) + c,,{s,t)ReFi''''''-'\t) 
+c,,{s,t)ReF^'''"''-'\t) , (16) 
Us,t)ReFi'''''''\t) + Us,t)ReF}'^''''-'\t) 
+c,4(s,t)ReFf ■^'^■-'^(t) + Q3(s,t)ReFf '^'-"^t) 
+c,,{s,t)ReFi''''''-'\t) . (17) 



The functions Cjj are polynomials of the Mandelstam variables; their explicit expressions 
are collected in Appendix El while the Laurent coefficients of the form factors F^"^^'^^^ 
{i = 1,2) were introduced in Eq. (fT^ . Even if the form factors are real for a physical 
(space-like) t, we write ReFl^'^^'^^\t) {i = 1,2) in the equations above for convenience of 
later use. 

The contribution to the differential cross section of the interference of the diagram in 
Fig. |21-(c) with the tree-level amplitude can be obtained from the contribution of diagram 
I21-(a) by replacing p2 — ps. This is equivalent to exchange s ^ t, so that one finds 



^^(y,ph Vertices) 



(phIrr.Ver.,2-c) 



a 
s 



1 

7t 



(2i,Irr. Ver.) 



{t,s) 



18) 



According to the definitions in Eqs. p5|l()P17j] . the r. h. s. of Eq. fll8|l involves the 
functions F'f'^'^^\s), which develop then an imaginary part above threshold (s > Am?). 
These imaginary parts do not contribute to the differential cross section at this order. 

Finally, the total contribution of the four diagrams in Fig. |21to the Bhabha scattering 
differential cross section at order is 



da. 



(V,ph Vertices) 



dn 



(phlrr.Ver.) 



2^ 

s 



(2i,Irr. Ver.) / ,n T^(2«,Irr. Ver.) 



[s,t) + V^- 



(19) 



The first, second, and third term within squared brackets in the r. h. s. of Eq. (fTI?|) are 
the s-s, t-t, and s-t channel interference amplitudes, respectively. 



4 Reducible Two-Loop Vertex Corrections 

In this section, we consider the interference of the two-loop reducible diagrams, shown 
in Fig. ini with the tree-level Bhabha scattering amplitude. The two reducible diagrams 
contain a one-loop vertex correction in both fermionic currents. Their contribution to 
the differential cross section at order can be written in terms of the one-loop UV- 
renormalized vertex form factors -F/^''* {i = 1,2) (see I26j) and of the functions Cij 
{i = l,2,j = 1,6) introduced in the previous section. It must also be observed that the 
term linear in {D — 4) is needed in the Laurent expansion of the one-loop vertex form 
factors. 
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(a) (b) 

Figure 3: Two-loop reducible vertex diagrams contributing to Bhabha scattering. 



The contribution of diagram 01- (a) to the cross section is given by 



da. 



(V,ph Vertices) 



dn 



a 



(phRed. Ver.|21-{a)) S 



i^^(2^,Red.Ver.)^^^^) + 1 ^-•) (s, t) 

st ij 



(20) 



where the Laurent expansion of ^(2Z,Red. ver.) 



^^(2/,Red. Ver.,-2) ^^(2Z,Red. Ver.,-1) 

(2«,Red. Ver.) 



with 



(2i,Red. Ver.,-2) 



(D - Ay 



{s,t) = Cii{s,t) 



(2«,Red. Ver. -1) 



{2Z,Red. Ver.,0) 



(D-A) 



(2«,Red. Ver.,0) 



+ 0{D-A), (21) 



{S,t) = Ci2{s,t) 



ReF("--^)(t))'-(lmF("'-^)(t) 

+2Qi(s,t)[ReFi('''"'^(t)ReFf''°^(t)-ImFf''"'^(t)ImFf''°^(t) 
+2q3(s, t) [ReFf (t)ReFf ''^^ (t) - lmFi^''-^'> {t)lmF^^''^^ (t) 



[22) 



(23) 



(S,t) = Q4(s,t) 



(ReF(^''-^)(t))'-(lmF(^''-^)(t) 

+2q2(s, t) [RcFI"'"'^ (t)ReFf '-^^ (t) -ImFi^'''"'^ (t)ImFi^"'°) (t) 
+2ci6(s, t) [ReFf ''"')(t)ReFf ''°)(t) -IniFf ''"'^(t)ImFf ''°)(t) 
+2ca(s, t) [ReFi('''"'^(t)ReFf '''^(t) -ImFi('''"'^(t)ImFi('''')(t) 
+2q3(s, t) [ReFf ''"'^(t)ReFi''''^(t) -ImFf '■"'^(t)ImFi''''^(t) 



+Cil(s,t) 



fReFf'-°)(t))'-flmFf''°)(t) 



2 / „^ N 2 

+2ci3(s, t) ReFf (t)ReFi'''°) (t) - IniFi^'''^^ (t)ImFf ''^^ (t) 



+Cj5(s,t) 



(ReFf'°)(t))'-(lmFf'°)(t) 



(24) 



Similarly, the contribution of diagram El- (b) to the cross section at order is given 



by 



^^(V,ph Vertices) 



(ph Red. Ver. 121- (b)) 



a 

S 



l^^(2^,Red.Ver.)^^^^^ + ^^(^'.Hed. Ver.) 



• (25) 



It is then possible to conclude that the interference of the diagrams in Fig. (21 with the 
tree-level amplitude gives the following contribution to the cross sections 



^^{V,ph Vertices) 



dn 



(ph Red. Ver.) 



a 

S 



1 ^^(2^,Red.Ver.)^^^^^ + l^f ^-■) t) 



(2Z,Red. Ver.) 



>,t) + vl 



(2Z,Red. Ver.) 



■ (26) 



The first, second, and third term within squared brackets in the r. h. s. of Eq. (j^Uj) are 
the s-s, t-t, and s-t channel interference amplitudes, respectively. 



5 Interference of One-Loop Vertex Diagrams 

In this section, we obtain the Bhabha scattering cross section at order deriving from 
the interference of the diagrams in Fig. ^ 

We begin by considering the contribution of the amplitude of diagram El-(a) squared; 
one finds that 

1 

m =—^h(s,t). (27) 

(phVer.Ver.|51-(a)) S 

In Eq. (j27j) . we introduce the function Ji that has the following Laurent expansion in 
powers of {D — 4): 



^^{V,ph Vertices) 



dn 



h{s,t) 



(D-4)2 (D-A) 



with 



1 



1 



C2l(s,t) 



2^22(5, t) 



2 / „, ,N N 2 

2 / ^1, 1^ \ 2 



(ReFi''^-'\t)) + (lmF(^''-^)(t) 



+C2i(s,t) ReFi'''-'>{t)ReFi''''"{t) + lmF^"'-'> {t)lmFi^''^\t) 



+C23(s,t) ReFi^''~^\t)ReF^''''"{t) + lmF'i'''''>{t)lmF^''''"{t) , (30) 



,(1«,0) 



/f)(s, t) = C22{s, t) \^eFl''^-'>{t)ReFl'''''>{t) + ImFi^"'-^^(t)ImFi^"'^^(t) 
+C2i(s,t) [ReFf '■~')(t)ReFf '■'^(t) + ImF^ ''"')(t)ImFi('''')(t) 
+C22(s,t) [ReFf ''"')(t)ReFf '■°)(t) + ImFf ''"')(t)ImFi'''°)(t) 
+C23(s,t) [ReFf '■"'^(t)ReFi''''^(t) + ImF^ ''"'^(t)ImFi'''')(t) 



.(i«-i) 



(l/,0). 



.(1«,0), 



(2^ 



(29) 
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(a) (b) (c) (d) 

Figure 4: One-loop vertex diagrams contributing to Bhabha scattering. 



(ReFi'^''\t)Y + (lmFi'^''\t) 



+ ^C2l{s,t) 



+C27is,t) 



(31) 



The expression of the functions Cij in terms of the Mandelstam invariants can be found 
in Appendix El while the expression of the one-loop vertex form factors appearing in the 
functions Ji can be found in ^12,2^ and are collected in [26j. 

Also the contribution of the square of the diagram in Fig. |3J-(c) can be easily expressed 
by employing the function Ji: 



da. 



(y,ph Vertices) 



dn 



(phVer.Ver.0]-(c)) 



1 



s 



hit, s). 



(32) 



The interference between the diagrams in Fig. 01- (a) and Fig. El-(b) generates a term 
in the cross section that can be written as 



^^(y,ph Vertices) 



dn 

where the function I2 is given by 



- ^1 

(phVer.Ver.El-(a)-(b)) S t"^ 



h{s,t) 



h{s,t) 



r(-2), 



S,t) 



+ 



(D-4)2 (D-A) 



+/f)(,,t) + 0(D-4) 



with 



4-) 


> 


t) 




It' 


> 


t) 


= 2lt'\ 


If) 


> 


t) 


= 2l['\s 



X 



ReF. 



(i«,o) 



it) 



ImF. 



(i«,o) 



(t) 



(33) 

(34) 

(35) 
(36) 

(37) 



Similarly, the interference between the diagrams in Fig. 111(c) and Fig. 111(d) generates 
the following contribution to the differential cross section 



da. 



(Viph Vertices) 



dn 



m = Kl2{t,S). 

(phVer.Ver.El-(c)-(d)) S 



(St 
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The interference of the diagrams in Fig. 01- (a) and Fig.01-(c) can be expressed in terms 
of a third function, J3: 



^^{V,ph Vertices) 



a 1 ^ , 
m = his,t). 

(phVer.Ver.gl-(a)-(c)) S St 



(39) 



Also in this case, it is convenient to exphcitly write the Laurent expansion of the function 



r(-2) 



{S,t) 



r(-l), 



+ 



s,t) 



with 



(D-4)2 {D-4) 
Cn{t,s)Fi'''-'\t)ReFi''^-'\ 



(40) 



(41) 



s,t) = cu{t,s)Fi^''-'\t)ReFi 



+cii(t, s) Fi^"'~'^(t)ReFf '•°^(s) + F^'"'"' {t)ReF^ 



,(1«,0), 



+ci3(t, s) Fi'''-'\t)ReF^'''">{s) + F^'"'"" {t)ReF^ 



{11,0), 



^{11,0), 



s) 



(42) 



/f(s,t) = c,,{t,s)Fr'-"meFi 



.(1^-1)/ 



+ci2(t, s) Fl'''~'\t)ReF^'''">{s) + F^'''"' {t)ReF^ 



{11,0), 



.{11,0), 



.{11,-1), 



+Cll(t, s) 



Fi'''-'\t)ReFr'>{s) + Fr'>{t)ReF^' ^ ^(s) 



+ci6(t, s)Fi('''-'^(t)ReF2^"'^^(s) + ci6(s,t)F2^"'^^(t)ReFi^ 



.(1«,0) 



{11,0), 



Ci3(t,s)F(^''~'Ht)ReF2^^''^^(s) + Ci3(s,t)F2^^''^^(t)ReF,^ 



ci3(s,t)Fi^''°)(t)ReFr''^^(s) + ci3(t,s)Fr''^^(t)ReF, 



.(i«,o) 



(li.O), 



.(i;,o). 



+cir(t, .)Fi^''°^(t)ReF2^^''^^(.) + cn(t, .)Fr''^^(t)ReFi^ 



(1«,0), 



.(li,0). 



(43) 



Finally, it is easy to verify that all of the interferences between s- and t-channel 
diagrams in Fig. lUgive a contribution to the cross section identical to the one in Eq. (jSHI), 
while the squared amplitude of the diagram lll-(b) (jll-(d)) coincides with the squared 
amplitude of diagram (a) (Ill-(c)). We can then conclude that the total contribution of 
the interferences of the diagrams in Fig. |3]to the Bhabha scattering cross section is given 
by 



^^(V,ph Vertices) 



dQ 



(phVer.Ver.) 



a 

S 



St 



+ -{2h{t,s) + l2{t,s)) 



(44) 



6 Interference of One-Loop Vertex and One-Loop Box 
Diagrams 

In the present section, we consider the interference of the diagrams in Fig. 0] with the 
diagrams in Fig. 
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(a) 



(b) 



(c) 



(d) 



Figure 5: One-loop box diagrams contributing to Bhabha scattering. 



In order to write down the contributions that these interferences provide to the 
Bhabha scattering cross sections at order a^, it is convenient to introduce the functions 
J4^j(xi, X2, X3), where i = 1, . . . ,3 and where Xj (j = 1, . . . , 3) represents one of the Man- 
delstam invariants s,t, and u; these functions have the following Laurent expansion in 
(D-4): 



hAxi,X2,X3) 



lii'^\xi,X2,X3) lii^\xi,X2,X3) 



+ /f/(a;i,a;2,X3) + 0(D-4). (45) 



{D - 4)2 (D-A) 
The coefficients of the Laurent expansion are given by 



Ili '{Xi,X2,X3} 



li/\xi,X2,X3) 



Fi''-'\x,)) Br-'\x2,X3) 



xsRe 
XsRe 



(46) 



f[''-'\x,)X Br^'\x2.X3) + (Fr\xi)Y Bi'^-'\x2,X3) 



(47) 



4?(a;i,a;2,X3) 



X3Re^[Fi'^-'\x,)y Bi'^''\x2,X3) + (^^^^1))* ^1"'°^ 
+ {Fr\-i)yBi^^-^\x2,xs) + (Ff ^°)(xO)*5fff (-.,-3) 



2^2, 3:3) 



(4J 



Besides for the one-loop vertex form factors already employed in the previous sections, 



one encounters in the equations above the functions B, 



1,...,3, J = -1,0,1) 



introduced in J2] to describe the contribution of the diagrams in Fig. El to the Bhabha 



scattering cross section^ at 0{a^). The auxiliary functions B^ 



4,...,6,j 



-1,0) 



are introduced for the first time here, and their expressions in terms of HPLs and the 
dimensionless variables x, y, and z can be found in pEj . 

The interference of diagram^ (a) with diagram El- (a) gives the following contribution 

^While the analytic expressions of the functions and iJ^^^''"-* are explicitly given in the 

expressions of the functions are not. They are the coefhcients of {D — 4) in the Laurent expansion 

of the functions i?!^'"* (Eq. (45) in ^2]), that were not needed in that context. Their expression in terms 
of HPLs and dimensionless variables can be found in • 
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to the differential cross section: 

^^(Vjph Vertices) 



— ^ I { 

(phVer.Boxlll-(a)-El-(a)) ~ Ist^ ^'^'^ ' ' 



(49) 



Tlie notation lias been chosen to clarify that Eq. ()49|) is related to the interference of two 
t-channel diagrams. 

Similarly, it is possible to write all the interferences of diagrams 01- (a) and0]-(c) with 
the diagrams in Fig. in terms of the functions I^^i. 



7 {V,ph Vertices) 
(10 2 




0? \ ^ , , 

'-±b i 


(50) 


Ui i 


(pn Ver.Boxl^h- (a) -|OJ- (b) ) 


^^(Vjph Vertices) 




1 


(51) 




(ph Ver .BoxEl- (a) -El- (c) ) 


^^(Vjph Vertices) 






(52) 




(ph Ver .BoxEl- (a) -0- (d) ) 


^^{V,ph Vertices) 




1 ^ , , 

^-/4,l(.,.,t), 


(53) 


dn 


(ph Ver .BoxEl- (c) -El- (a) ) 


^^{V,ph Vertices) 




\ ^ , , 

^-/4,3(3,«,t), 


(54) 


dn 


(ph Ver .BoxjH- (c) -0- (b) ) 


^^{V,ph Vertices) 




1 

— ^/4,2(s,t,s), 


(55) 


dn 


(ph Ver .BoxjH- (c) -El- (c) ) 


^^{V,ph Vertices) 




0? 1 


(56) 


dn 


(phVer.Boxlll-(c)-El-(c)) 



The interference of diagram El-(b) (01-(d)) with a diagram in Fig. El is identical to the 
interference of diagram 01- (a) (01- (c)) with the same diagram in Fig. El Therefore, the 
total contribution of the interferences between one-loop vertex diagrams and one-loop 
box diagrams to the Bhabha scattering cross section is twice the sum of Eqs. 
namely 



^^(Vjph Vertices) 



dn 



(ph Ver. Box) 



a 

{.h,l{t,t,s) +l4^s{t,U,s) +l4^i{s,S,t) +l4^3{s,U,t)) 

st 



+^ {h,2{s, t, s) - 14,2(5, u, s)) 



(57) 



7 Interference of One-Loop Box Diagrams 

With reference to Eq. ()10|) . it is possible to further split g-(^'P*^^°^®^) the sum of two 
contributions; the first originates from the interference of two-loop photonic box diagrams 
and tree-level diagrams, the second originates from the interference of two one-loop box 
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diagrams: 



da. 



(y.ph Boxes) 



da. 



(V,ph Boxes) 



dn 



dn 



+ 



^^(y,ph Boxes) 



(ph Box Box) 



dn 



(58) 



(2LBox) 



The two-loop photonic box diagrams in the m 7^ case are at the moment still unknown, 
and the calculation of the second term in the r. h. s. of Eq. (j3Hj) currently remains as an 
open problem (see [201 HI])- 

On the contrary, the calculation of the first term in Eq. (j58p is, in principle, straight- 
forward. The interference of every pair of the one-loop box diagrams shown in Fig. El 
provides a contribution to the differential cross section of the form 



da. 



(V,ph Boxes) 



dn 



a 



(ph Box BoXjij) 



—Cij{x,y,z) 



(59) 



where the indices i,j run over the diagram labels = a,b,c,d) and where we have 
introduced new functions Cij. The explicit expression of the latter, already continued to 
the physical region s > 4m^, is particularly long and can be found in 



8 Soft-Photon Emission at Order cx^ 

All of the two-loop photonic corrections discussed in the previous sections are UV renor- 
malized, but they still include double and single poles in {D — 4). These singularities 
have an IR origin, and they can be eliminated by adding the contribution of the real 
soft-photon emission diagrams at order to the virtual cross section. 

Before discussing the soft corrections to the Bhabha scattering differential cross section 
of order a^, the reader is reminded that the soft corrections at order a^, discussed in detail 
in can be written in the factorized form 

/a\ daf{s,t,m'^) / a\ da^{s,t,rn?) 

— dn dn ' 

where o"^(s, t, m^) is the tree-level cross section obtained by calculating the traces of Dirac 
matrices in D dimensions and where 5*1^ is defined as 

4 

-S'lR = 4 Jij , Jij = ej {pi ■ Pj) hj , (61) 
i=i 

with ei = €4 = 1 and €2 = €3 = —1, and 

1 m(^-4) r d^~^k 1 

" r (3 - f ) 7r(^-4)/2 47r2 J ko {pi ■ k) {p, ■ k) ' ^ ' 

The integral in Eq. (jU^ can be found in (see also Appendix A of 13 ); the integration 
over the momentum of the soft photon {k) is restricted to the region \k\ = k^ < uj, where 



13 








(a) 



(b) 



(d) 



(e) 



(f) 



(g) 



(h) 



Figure 6: Examples of double photon emission from tree-level diagrams and single photon 
emission from one-loop diagrams. 



uj is the cut-off on the energy of the unobserved soft photon. The expansion of (Jq{s, t, m?) 
in powers of {D — 4) is 



daQ{s,t,m'^) 

dn 



dn 

+ p-4)2 



daj^\s, t, m?) 

dn 



dn 

+ 0{{D-Af) 



(63) 



where cro(s,t,m^) is the well known tree-level cross section, and 



dal^\s,t,m'^) o? 

dQ s 

daf\s,t,m?) 1 

dVt s st 



+ - 



+ ■ 



st 



\{s+tf-\st-m\s+t) 



:St 



(64) 
(65) 



The contribution of the s- and t- channel diagrams, and of their interference to a^^^ (s, t, m?) 
and al^\s,t,m'^), is evident in Eqs. ()64I65|) . 

There are two different kinds of soft-photon emission diagrams contributing to the real 
corrections to the cross section at order a^: 

i) the tree-level diagrams with the emission of two soft photons (some diagrams be- 
longing to this class are shown in Fig. |Hl-(a)-(d)), and 

ii) the diagrams which include a one-loop correction and the emission of a soft photon 
from one of the external legs (see the examples in Fig. ini-(e)-(h)). 

Since the soft-photon corrections in QED exponentiate, the contribution of the double 
emission diagrams to the Bhabha scattering differential cross section is given by 



da2^'^°'^^^^\s,t,m'^) 1 /a\^ daQ{s,t,m'^) 2 
n) dn " 2 W J dn ™ 



(66) 



The emission of a single photon from one-loop diagrams, interfered with the tree-level 
single photon emission graphs, amounts to 



a\2 daf'''"^''^(s,t,m2) /ax 2 da^^'''\s,t,m^) 



TT 



dn 



TT 



dn 



IR : 



(67) 
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where crj^'^^ is the UV-renormahzed virtual cross section at order a^. The superscript D 
was introduced as a reminder that the Laurent expansion of a^'^^ must be known up to 
and including terms linear in {D — A); the reason for this is that these linear terms give 
rise to a finite contribution when multiplied by the single pole present in Si^. The real 
corrections of order which originate from single-photon emission diagrams and which 
include a fermionic loop (as, for example, diagram |B]-(h)), were calculated in [13J; they 
include single IR poles that cancel against the virtual corrections of order that also 
include a photon self-energy insertion. This set does not play any role to the present 
discussion and is thus systematically ignored. 

With the above in mind, the one-loop virtual cross section appearing in Eq. (jU7j) can 
be written as 



dQ 



dn 



+ 



da^^'^\s, t, im?) 



(iiy) 



dn 



(6^ 



where the subscript V indicates the contribution of the interference between vertex graphs 
and tree-level amplitude; B stands for the cross section generated by the interference 
between one-loop boxes and tree-level diagrams. Furthermore, we split the terms on the 
r. h. s. of Eq. ()68p as follows: 



da'i''^\s, t, im?) 



dn 



daY {s, t, rn?) 



dn 

+0{iD 



da\ 



dn 



(69) 



with j = V,B. The singular and finite parts of the one- loop virtual cross section, corre- 
sponding to the first term in the r. h. s. of Eq. (jUHj) . are well known. Their expression in 
terms of Mandelstam invariants and HPLs can be found in Eq. (43) and Eq. (49) of jT2] • 
The terms proportional to {D — 4) which arise from box and vertex one-loop corrections 
are given by 



d(7^^'^\s, t, m^) 



dn 



{ll,B) 



a 
s 



ReBi''''^ («, s)) + ^ (Re5f -'^ (s, t) + ReB[''''^ (t, s) 



-5^''''^(n,t) + Re5f'')(u,s; 



(70) 



and 

da'l^'^\s, t, m^) 



dn 



(iiy) 

respectively, with 

V}''^'\s,t) = Ca(s,t)ReF,('^''^(t)+Q2(s,t)ReFi('''°Ht)+ Q4(s,t)ReFf''-')(t) 



[71) 



+Q6(s,t)ReFi'''°)(t) + Q3(s,t)ReF2^"'^^(t) 



(z = l,2). (72) 
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All of the functions appearing in the r. h. s. of Eqs. ()70ll72j) were introduced in the previous 
sections, and their expressions in terms of Mandelstam invariants and HPLs have been 
collected in 

As is mentioned above, the calculation of the integrals Iij which appears in Eq. (jF)T|) 
has been carried out in [T^, up to terms linear in {D — 4) excluded. At first glance, 
it appears that the calculation of such terms is needed, since they provide, in the limit 
D —>■ A, a non vanishing contribution to both Eq. (jUBj) and Eq. ()67|) . However, it is 
possible to prove that this is not the case. In order to proceed with our proof, one needs 
to split the one-loop UV-renormalized virtual corrections in an IR-divergent part and a 
finite reminder, as is exemplified in the following: 



daY {s, t, rn? 



da[ 



dn 



The Laurent expansion of Sir has the form 



+ 



da[ 



{n,j) 



dn 



(73) 



S- 



IR 



(-1) 

IR 



si^ + {D-4)S[^ + 0{iD-Ar) . 



:(i) 



(D-A) 



(74) 



The cancellation of IR divergencies in the order cross section guarantees that" 



da 



(v,-i), 



s, t, m 



dn 



dal 



s, t, m 



(ll,B) 



dn 



dao{s,t,m'^) 

+ ~ >Jt" 



(iiy) 



dn 



0. 



(75) 



By employing Eqs. ()66l67p in combination with Eqs. ()63I73I74|) . one can prove that the 
non-vanishing term proportional to S^^ appearing in the double emission cross section 
(Eq. (jnSl)) is 



^^(S,double) ^g^f^^2^ ^ d(To{s, t,m'^) (_i) (1) 



dn 



dn 



(76) 



while the non vanishing term proportional to appearing in the single-photon emission 
cross sections at order is 



da. 



(5,slngle) 



dn 



da^' ^\s, 1,171^ 



dn 



+ 



da^^'-^\s,t,m') 



dn 



{iiy) 



S- 



(1) 

IR 



(77) 



Therefore, we can conclude that, due to Eq. (ffH|) . the non vanishing terms proportional 
to cancel out in the total real emission cross section at order a^, given by the sum of 
Eqs. (IMIHTjl . 



9 Cancellation of the IR Singularities 

The IR divergencies in the real corrections at order (Eqs. (j66l67|) ) should cancel the 
IR singularities present in the virtual corrections discussed in Sections OHZI and the ones 
arising from the interference of the (yet unknown) two-loop photonic box graphs. 

■^We remind the reader that the photon self-energy diagrams are IR- finite (see |13p. 
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Figure 7: Example of the cancellation of the IR divergencies in the two-loop irreducible 
vertex corrections (first line) and in the interference of the one-loop direct box in the 
s-channel with itself (second line). 




Figure 8: Residual IR poles proportional to (^(2) in the s-channel cross section. 



As in the case of the cancellation of the IR divergencies of the virtual cross section at 
order a^[Np = 1) discussed in it is possible to organize the contributions to the cross 
section in IR-finite blocks by pairing the virtual corrections originating from a certain set 
of diagrams with an appropriate subset of the soft-photon emission corrections. As an 
example, in the first line of Fig. [7| we illustrate the cancellation of the IR poles present 
in the interference of the graph in Fig. |21-(a) with the t-channel tree- level diagram; in all 
the terms in the 1. h. s., the product of two graphs represents the contribution of their 
interference to the Bhabha scattering differential cross section. Once again, the gray circle 
represents the sum of the UV-renormalized two-loop photonic vertex corrections to the 
electron current in the t-channel photon-exchange diagram. In the second line of Fig. [71 
we provide another example of the cancellation of the IR divergencies in the box-by-box 
sector. Similar relations can be found for all of the contributions to the cross section at 
order a'^. 

A special case is represented by the virtual corrections in which the coefficients of the 
IR poles, when calculated in the non-physical region s < 0, include HPLs with two or 
more zeroes in the rightmost positions in the weight list. When performing the analytic 
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Figure 9: Residual IR poles proportional to ({2) in the s-t-channel interference contribu- 

, ■ , (V.ph Vertices) 

tion to 02 



continuation to the physical region s > 4m^, these HPLs generate real terms proportional 
to C(2), where ( is the Riemann Zeta function. For example, by replacing the non physical 
dimensionless variable x according to x — > —x' + ie, one finds that 

H{0, 0; x) H{0, 0; -x' + ze) = H{0, 0; x') - 3C(2) + iTTH{0; x') , (78) 

or 

H{-1, 0, 0; x) H{-1, 0, 0; -x' + ze) = -H{1, 0, 0; x') + 3C(2) - i7iH{l, 0; x') . (79) 

The part of the IR pole, that is proportional to the C(2) factor arising from analytical 
continuation does not cancel in the combination with the soft-radiation contributions that 
eliminate the other IR singularities. This kind of behavior, already observed in p!3| in the 
discussion of the 0{a'^{Np = 1)) cross section, is again encountered in the cases illustrated 
in Figs. EHm 

In the total differential cross section, the residual poles of Fig. |S1 cancel themselves 
out. Combining the various contributions shown in Fig. |H1 the residual IR poles do not 
cancel. Residual single poles proportional to C,{2) also arise in the combinations shown in 
Figs. Uni and [TT] Clearly, such residual poles cancel once the contribution of the two- loop 
box graphs is included in the cross section at order a'^. The five two-loop photonic box 
topologies are shown in Fig. 

10 Logarithmic Expansions 

In order to check our calculations against the results available in the literature, the contri- 
butions to the cross section at order described in the previous sections are expanded in 
the limit in which the mass of the electron is small with respect to all of the Mandelstam 
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4(JnJi3 + Ji2^: 



13) 





2J| 



13 





+ ("^11 + ^12) ) ( X 




single IR 
pole oc C(2) 








-(^11 + ^12) ) ( X 




single IR 
pole oc C(2) 



Figure 10: Residual IR single poles proportional to ({2) in the interference of one-loop 
vertex by one-loop box diagrams. 




Figure 11: Residual IR poles proportional to C(2) in the interference amongst one-loop 
box diagrams. 
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(a) 



id) 



(d) 



id) 



Figure 12: Two-loop photonic box topologies. 



invariants s, t, and u^. In this limit, it is customary to write the two-loop photonic cross 
section as follows (see P HHj): 



da. 



(pt) 



dan 



dan 



4 

m^/ \ s 



(80) 



In the equation above, 
the previous sections. 



ph Vertices, ph Boxes. Following the notation adopted in 



da. 



(y,ph Vertices) 



da. 



{V,ph Vertices) 



dVt 



dn 



+ 



da. 



{V,ph Vertices) 



(ph Irr. Ver.) 



dn 



(phRed. Ver.) 



+ 



da. 



(y,ph Vertices) 



dn 



+ 



da. 



(y,ph Vertices) 



(ph Ver. Ver.) 



dn 



(81) 



(phVer. Box) 



while a^'"^^ Boxes) introduced in Section [7| The corresponding cTg'^'*'' cross sections 
were obtained by pairing virtual and soft-photon emission contributions as described in 
Sectionini In several points of our discussion, we stressed the fact that a non-approximated 
calculation of the contribution of the two-loop photonic box diagrams to the cross section 
is still missing. However, all of the coefficients (z = 2, 1, 0) in the expansion in Eq. (jHUj) 
are completely known; the first two can be found in [H], while 6^^ was recently obtained 
in ^Uj. Therefore, by employing the m^/s — limit of the results presented here in 
combination with jHl 02] , it is possible to indirectly obtain the m^/s limit of the 
contribution of the yet unknown two-loop photonic boxes (and corresponding soft-photon 
emission corrections). In Appendix ^ we report the expression of such a contribution, 
both before and after adding the corresponding soft-photon corrections. Moreover, the 
expansions in the m^/s limit of all the contributions to the Bhabha scattering cross 
section discussed in the present paper can be found in [22] • 

It is known that the small-angle Bhabha scattering cross section is completely deter- 
mined by the Dirac vertex form factor [4J. In particular, one finds that for the virtual 
cross section 



dan 



82 



where F^^^''^ and are the UV-renormalized vertex form factors already employed in 

this paper. The IR poles present in the form factors are easily removed by adding the 

''Note that this expansion is not vahd for very smaU scattering angle, corresponding to |<| < m^, and 
for almost-backward scattering, corresponding to |u| < m^. 
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soft emission contributions. By introducing the IR-finite form factors 



Fi''\t) = Frit) + Jn + Ji3, 



.{11) I 



(83) 



one finds that 



da. 



(Ph) 



0^0 



dan 



{2l,ph) 



;i-^+e 



2^2 



In^ 



93 . „ 2 ^4cj2 



6C(3)-3C(2)-- + 91n(0-41n 



[1 - In (0] 



- 21n( — 1 [7-61n(0] 



- 9C (3) + -C (2) - 12C (2) In (2) - (2) 
4 5 



27 93 9 ^ 

+ y + 6C (3) In (0 - 3C (2) In (0 " f 1^ (0 + | In iO" 

+ InM — ) [2-41n(0 + 21n2(0] 



+ In 



s 



[8-141n(0 + 61n2 (^)] 



The variable ^ is defined as 



cos 6* 



^4) 



^5) 



with 6 the scattering angle in the c. m. frame. By expanding the exact photonic correc- 
tions involving vertex diagrams in the 6^0 limit and by neglecting terms proportional 
to the electron mass, we recover the expression in Eq. (jH3), which agrees with ^U] and 
the theorem in jl] . Consequently, this represents a non-trivial test of our calculation. As 
expected, we observe that the interference between one-loop box diagrams and one-loop 
vertex corrections does not contribute to the small-angle cross section. The interference 
of one- loop box diagrams amongst themselves has a non-zero small angle limit, and all 
of the residual terms are proportional to C(2)- These residual terms cancel out once the 
contribution of the two-loop photonic boxes is added; for this reason they are excluded 
from the present discussion. The agreement with the results of PU] is clarified in Fig. IT^ 
where we plot as a function of the scattering angle 6 the result of ^D] and the corrections 
of order originating from vertex graphs (Sections EHHl plus corresponding soft emission 
contributions). It is easily seen that, at small angles, the vertex corrections completely 
determine the cross section. 

In addition, the expansion of the interference of one-loop boxes provides another strong 
test of both our calculation and the one discussed in JUl • The interference of some pairs 
of the diagrams in Fig. El gives origin, in the ratio da2/dao, to terms proportional to 



ln"e 



n 



m 



1,2 



(86) 



21 




Figure 13: Order Bhabha scattering differential cross section divided by the cross section 
in Born approximation, as a function of the scattering angle 9. The continuous line 
represents the result of ^101, while the dashed line represents the corrections involving at 
least a vertex graph. The beam energy is chosen equal to 0.5 GeV and the soft-photon 
energy cut-off ui is set equal to E. 



It is possible to observe that, in the sum of all the one-loop box interferences, such terms 
cancel out, and that they do not appear in the complete photonic cross section at order 

The photonic corrections to the Bhabha scattering cross section are now known up to 
terms of 0{rn? /s) excluded; it is possible to use the corrections calculated exactly in this 
work to estimate the relevance of the 0{m? / s) terms. We define 

dVL ~ dVL ~ dVL 

where the index i = Vertices, Box Box represents the contributions discussed in Sec- 
tions OHEl and respectively. In Figs. [HI and [13 we plot, as a function of the beam 
energy, the quantities 

It can be seen from these plots that the terms proportional to the electron mass become 
negligible for values of the beam energy that are very small with respect to the ones 
encountered in practically all of the e^e~ experiments. It is also reasonable to expect 
that the terms proportional to the electron mass are negligible in the corrections due to 
the two-loop photonic boxes. In this sense, the approximated cross section obtained in 
jni [ini [121 [in] should be sufficient for all phenomenological studies. 
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Figure 14: -Dvertices (is a function of the beam energy, for 9 = 10° (solid line) and 9 = 90° 
(dashed line). The soft-photon energy cut- off is set equal to E. 
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Figure 15: Dbox box OjS a function of the energy, for 9 = 10° (solid line) and 9 
(dashed line). The soft-photon energy cut-off is set equal to E. 
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Figure 16: Comparison of the vertex corrections with the complete photonic cross section 
at order a^. E and uj as in Fig. ■ 



Finally, in Figs. El and El it is possible to compare the various contributions to the 
Bhabha scattering differential cross section known at present with the complete photonic 
cross section in the m? / s limit [K]. The dashed line in Fig. ITHl corresponds to 
the contribution to the cross section defined in Eq. (|87p for i = Vertices, plotted as a 
function of the scattering angle. The local minimum at 6' ~ 80° and the maximum in the 
backward direction are due to spurious terms proportional to the monomials 

C(2)W(4), C(2)ln(^)lne, and C(2) In (4) - • (89) 

These terms are not present in the complete cross section (they cancel out against anal- 
ogous contributions deriving from the interferences of two-loop box diagrams with the 
tree-level amplitude and from the interference one-loop box diagrams amongst them- 
selves). Removing them from the vertex corrections, one obtains the dashed-dotted curve 
in Fig. El which is smoother than the dashed one. The solid curve represents the complete 
photonic cross section. As is already observed above, the vertex contribution reproduces 
the full result in the small-angle region. 

The dashed line in Fig. El represents the quantity 

^>h Boxes) ^(ph) ^(ph vertices) Box Box) ^>h2LBox) 

= -T —i = ^ + -h • (90) 

auo cLctq aao cicxq aco 

In the equation above, (j(p'^^°^^°^) -^^g introduced in Eq. (|87p. while ^-(p'^^'"^"^) jg ^^j^e 
contribution to the cross section of the two-loop photonic boxes interfered with the tree- 
level amplitude. The expression of the latter is given in Appendix El Also in this case, 
the behavior at ^ ~ 80° and at 6' ~ 180° is dominated by the spurious terms of the kind 
of Eq. (jH^ : by removing them, we obtain the dashed-dotted curve. 
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Figure 17: Comparison of the box corrections with the complete photonic cross section at 
order . E and uo as in Fig. 



For completeness, in Fig.UHlwe plot the photonic [lOj and Np = 1 |13j contributions to 
the Bhabha scattering cross section at order a^. The dotted line represents the photonic 
corrections. The corrections of 0{a^{Np = 1)) (dashed line) have, for this choice of 
uj {uj = E), a.n opposite sign with respect to the photonic corrections. However, it is 
necessary to say that the 0{a^{NF = 1)) of include large contributions proportional 
to ln'^(s/m^) that cancel out once the contribution of the soft pair production is included. 
Of the latter, only the terms proportional to ln"'(s/m^) (n = 1,2,3) are known (see [H]); 
we checked that they cancel the ln^(s/m^) term of the 0{a^{NF = 1)) cross section. 
The dashed-dotted line represents the sum of the 0{a'^{Np = 1)) cross section with the 
known terms of the pair production corrections^. The solid line is the complete order 
QED Bhabha scattering cross section, including photonic, Np = 1 and pair production 
contributions. 



11 Conclusions 

In the present paper, we obtained analytic non approximated expressions for all of the 
photonic corrections to the QED Bhabha scattering differential cross section at second 
order {0{a^)), except for the ones deriving from the interference of two-loop photonic 
box diagrams with the tree-level amplitude; at present, the integrals necessary to the 

^The pair production, calculated in up to terms enhanced by In" (s/m^) included, depends upon a 
cut-off on the energy of the soft electron-positron pair. In the numerical evaluation of Fig. ^| we set 

ln(D) = iln(^l^) ' ^^^^ 
with ln(D) defined in |S] and f2 numerically equal to the soft-photon cut-off: — uj. 
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Figure 18: Photonic, Np = 1, and total contributions to the cross section at order a^. E 
and uj as in Fig. The pair production cut-off is set equal to uj. 



calculation of the latter are not known. The calculations were carried out by retaining 
the full dependence on the electron mass m. The results are valid for arbitrary values of 
the c. m. energy s and momentum transfer t. 

We included a discussion of the soft-photon emission at order a^, employing dimen- 
sional regularization to handle the IR-divergent terms. We proved that the term pro- 
portional to {D — 4) in the Laurent expansion of the photon phase-space integral in 
Eq. does not contribute to the Bhabha scattering cross section. 

After subtracting the IR singularities by adding the contribution of the soft-photon 
emission graphs, we expanded the contributions to the cross section discussed in Sec- 
tions OHZI in the m^/s — >■ limit. In this way, it was possible to cross-check large parts 
of the result of jT^, as well as to efficiently test our calculations. By subtracting the 
contributions to the cross section obtained in this paper from the cross section of jTU], it 
was also possible to indirectly obtain the contribution due to the interference of the two- 
loop photonic diagrams with the tree-level amplitude, up to terms suppressed by positive 
powers of the electron mass. 

By comparing the non approximated results with the corresponding m? / s —>■ limit, 
we explicitly checked that the contribution to the cross section of the terms proportional 
to positive powers of the ratio m'^/s is negligible at high- and intermediate-energy e~^e~ 
colliders. It is reasonable that the same conclusion applies to the contribution to the order 

cross section involving the two-loop photonic box diagrams. For what concerns phe- 
nomenological studies, the results of [6iEllHE] therefore provide a complete expression 
of the virtual and soft-photon emission corrections to the Bhabha scattering cross section 
at order a^. Nevertheless, a future calculation of the two-loop photonic box graphs for 
m 7^ would represent a very interesting result in the field of multi-loop calculations. 
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A Cij Functions 

In this appendix, we collect the explicit expressions of the Cij functions employed in the 
main text: 



cii(s,t 

Ci2(s,t 
Ci3(s,t 

cu{s,t 

Ci5(s, t 
Cl6(s,t 

cn{s,t 

C2l(s,t 

C22(s,t 

C23is,t 
C24(s, t 

C25('5, t 
C27(s,t 



{s + tf - Am' , 

- [{s + tf - st-2m'{s^t 
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st + 



st' 
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A 8 
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2 ' 

t 
2 ' 



st 



st' 



sH 



+ -i 
Am' Am' A 



(92) 
(93) 

(94) 
(95) 
(96) 
(97) 
(98) 

(99) 

(100) 

(101) 
(102) 

(103) 
(104) 
(105) 



B Photonic Double Boxes 

We provide here the expansion of the interference of the two-loop photonic boxes with 
the tree-level amplitude in the limit m' / s — >■ 0. Below one can find the contribution of 
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the two-loop photonic boxes to the virtual cross section as well as the same contribution 
after the subtraction of the corresponding soft-photon corrections. 



Without including soft radiation 
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Including soft radiation 
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